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Abstract:We discuss damping of the collective rotational motion in A ~ 150 superdeformed nuclei by means 
of a shell model combining the cranked Nilsson mean-filed and the surface-delta two-body residual force. It is 
shown that, because of the shell structure associated with the superdeformed mean-field, onset energy of the 
rotational damping becomes E x ^2 — 3 MeV above yrast line, which is much higher than in normal deformed 
nuclei. The mechanism of the shell structure effect is investigated through detailed analysis of level densities 
in superdeformed nuclei. It is predicted the onset of damping varies in different supedeformed nuclei along 
with variation in the single-particle structure at the Fermi surface. 
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1 Introduction 

In recent experiments, up to around 20 different rotational bands are observed in 
a normally deformed rare-earth nucleus, and several superdeformed rotational bands 
are identified in a single superdeformed nucleus M. The observed rotational bands 
usually lie in the region near the yrast line and are often interpreted in terms of in- 
dependent particle excitations in a rotating deformed mean-field. There are of course 
numerous excited levels at higher excitation energy, but these levels do not necessar- 
ily form rotational band structure characterized by sequences of strong E2 transitions 
since damping of the collective rotational motion is expected @. In the highly ex- 
cited region, say at 1 to 2 MeV above yrast line for normally deformed nuclei, the 
level density at a given spin is as large as 10 2 to 10 3 MeV -1 , and the corresponding 
level spacing becomes comparable with or smaller than the size of matrix elements 
(order of 10 keV) of the residual two-body interaction. The residual interaction then 
becomes effective to cause mixing of many-particle many-hole configurations in the 
rotating deformed mean-field potential. Since different configurations respond differ- 
ently to the Corioli force due to different single-particle alignments, the configuration 
mixing results in loss of collectivity or damping of the collective rotational motion. 
Experimental data on the damping of the rotational motion are accumulating for 
normally deformed rare-earth nuclei from the analysis of the quasi-continuum part 
of gamma-ray spectra containing gamma-rays emitted from the highly excited lev- 
els 0, f§. In particular, a fluctuation analysis of the quasi-continuum points to the 
presence of about 30 rotational bands in a nucleus, thus confirms occurrence of the 
rotational damping [[J. 

The mean-filed models which have been very successful for rotational bands near 
yrast line have to be extended in order to describe the rotational damping caused 
by the configuration mixing effect due to the residual interaction. Such extensions 
can be formulated by using the shell model diagonalization on the basis of many- 
particle many-hole configurations associated with the cranked mean-filed [5-7]. In 



fact, a realistic model which combines the cranked Nilsson single-particle basis and 
the surface-delta residual two-body force |F| gives results consistent with the data from 
the fluctuation analysis method for normally deformed rare-earth nuclei. It predicts 
that the rotational damping sets in at around 1 MeV above yrast line. 

The rotational damping in superdeformed nuclei is expected to differ significantly 
from that in normally deformed nuclei since the level density of superdeformed states 
is affected by the shell structure in the single-particle spectra^], HJ, which is the char- 
acteristic feature specific to the superdeformation. The shell-effect on the rotational 
damping is pointed out by an previous work |5[ for 152 Dy, which however adopted 
schematic residual interaction with constant matrix elements treated as parameters. 
Recently, the fluctuation analysis method has been applied to a superdeformed nu- 
cleus 143 Eu with a first decisive and quantitative data for the onset of the rotational 



damping in a superdeformed nucleus fllP] . It extracted an effective number of superde- 
formed rotational bands which is not very different from the effective number of bands 
in normally deformed nuclei. Since the shell structure in superdeformed nuclei is not 
uniform even in A ~ 150 region, systematic and quantitative analysis is required in 
order to reveal the characteristics of the rotational damping. 

The present paper analyzes systematically the rotational damping in superde- 
formed nuclei in the A ~ 150 region by means of the shell model diagonalization 
with realistic two-body residual interaction ||. In particular, the analysis is focused 
on the onset of the damping since it is most relevant to the data extracted from 
the fluctuation analysis method. In order to clarify the shell structure effect on the 
rotational damping, we also look into the level densities in detail in Section |3]. It 
is revealed that, as consequence of the shell structure effects, the onset of damping 
varies in different species of superdeformed nuclei. This is because the single-particle 
spectra as well as the shell-gap at the Fermi surface varies as N and Z change. 



2 Formulation 

The microscopic model adopted in the present study is essentially the same as 
the formulation for normally deformed nuclei |7j], except difference in the mean- field 
potential. In describing superdeformed states, suitable choice of the deformation 
parameters is important since the shell structure characteristic to superdeformation 
depends on the deformation. 

We start with the cranked Nilsson single-particle Hamiltonian 

h = /^Nilsson — Ujx (1) 

with u being the rotational frequency. The pairing potential is not introduced. For 
the Nilsson potential, the quadrupole and hexadecapole parameters e 2 and €4 are con- 



sidered as shape variables. We adopt the Nilsson parameters given by Ref . [O] , which 



gives better description of proton 1\ orbit than the parameter [E| used in Ref. 



The difference in the Nilsson parameters does not affect the following argument ex- 
cept detailed behavior. With reflection symmetry, the single-particle routhian orbits 
of the cranked Hamiltonian keep the parity and the signature quantum numbers. In 
order to guarantee smooth change in the single-particle spectrum as a function of the 
rotational frequency, an diabatic basis for the single-particle orbits is constructed in 
a similar way to Ref. JT3| |. 

An reference superdeformed configuration is defined for a given nucleus by putting 
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Table 1: The equilibrium deformation parameters 62 and £4 at / — 50h determined by Strutinsky 
method for each superdeformed nucleus. 



iV neutrons and Z protons up to the Fermi surface in the diabatic single-particle spec- 
trum. When the reference configuration is not unique due to crossings at the Fermi 
surface, we choose one by referring to the high-N configuration suggested from the ob- 



served yrast superdeformed band [14]]. Deformation parameters (62,64) are determined 
as functions of spin by minimizing the total energy of the reference configuration with 
the Strutinsky shell correction defined at a given spin I ||T5"|1 

E(I) = E micm (I) - E smooth (I) + E RLD (I). (2) 

Here the microscopically calculated total energy is given by 

£micro (/)= £ eiM+UjI (3) 
occupied i 

with the rotational frequency uij given 

Jx{ui) = JviM = 7 > ( 4 ) 

occupied i 

and Ci(uj) and j x i{oj) being the diabatic single-particle routhian energy and expectation 
value of j x . The Strutinsky smoothed total energy is given by 

E smooth {I) = e*(£) + ^ l (5) 

and 

Jx(v) = ^ixiiCb) = J, (6) 



with §i,j x i being e^j x i weighted by the smoothed occupation number hi |]16fl . 
The rotating liquid-drop energy E RLD (I) is calculated according to Ref. [17]. The 
deformation parameters thus determined are listed in Table. [I] at a representative 
spin / = 50h. The spin dependence of the deformation is weak in all the nuclei under 
consideration, e.g. (e 2 , e 4 ) = (0.577, 0.010) ~ (0.5700.030) for / = 10 ~ 70h in 152 Dy. 
An example of the calculated cranked Nilsson spectra with the optimized deformation 
is shown in Fig. [I] for 152 Dy. 




Figure 1: The single-particle routhian spectra for 152 Dy with the equilibrium deformation deter- 
mined by the energy minimization. The deformation parameters (€2,64) change continuously from 
(0.577,0.010) at u = 0.07(MeV/7i), I = 107i to (0.570,0.030) at w = 0.74(MeV/fr), I = 70h . 

In a superdeformed nucleus, there are usually more than two coexisting local minima 
of the total energy surface in the deformation parameter space, which correspond to 
superdeformed and weakly deformed mean- fields. If the mean-fields are stable (the 
minima are deep), excited states are expected to be built upon each mean-field. At 
high spins (I^AOh) where the superdeformed minimum becomes the lowest, states 
near the yrast line are dominated by excited superdeformed states. 

The basis for the excited superdefomed states is formed by many-particle many- 
hole (np-nh) configurations in the diabatic cranked Nilsson orbits which cover the 
energy range 5E = ±4.8MeV above and below the Fermi surface. The deformation 
parameters and rotational frequency specifying the single-particle basis for np-nh 
configurations are taken common to those for the reference configuration. Energy of 
the basis states is given by Eq. (Q), in which the microscopic energy is given for each 
np-nh configuration(labeled by /i) in the same way as in Ref. ||,z.e., 

^ic ro(/) = EM) + WjJ + _L_(j _ JxM) f (7) 

where 

E„{u I )= ]T ei(u)i) (8) 

occupied i 
in ft 



and 
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(9) 



is the dynamic moment of inertia of the configuration /i. The Strutinsky renormaliza- 
tion is important not only to determine the deformation but also to describe gamma- 
ray energies of the collective E2 transitions. The basis of the excited superdeformed 
states is denoted as with the energy {E^(I)}. It keeps good quantum num- 

bers with respect to the total signature and parity. Excited states built upon weakly 
deformed mean-fields are neglected in the present paper. At low spins (I^AOh), 
they may dominate the yrast structure and the superdeformed states are to be sur- 
rounded by these weakly deformed states, while coupling between the superdeformed 
and weakly deformed excited states is expected to be weak if there exists a potential 
barrier between the local minima. We neglect the presence of the weakly deformed 
states since the coupling is hardly treated by a shell model description. Hereafter 
the yrast line refers to an envelope of the lowest energy superdeformed states as a 
function of spin. 

For the residual two-body interaction, it is natural to assume that the effective 
nuclear force is not dependent on the deformation. We adopt the same surface-delta 
interaction (SDI) ]18[ with the same strength as used for describing the rotational 



damping in normally deformed nuclei so that results here can be compared with those 
for normally deformed nuclei 0. The strength of SDI is Vq = 27.5/A(MeV), which 



is taken from Ref. |19[]. The two-body interaction is decomposed into mean- field part 
associated with the reference configuration and its residual, and only the residual 
two-body force is taken into account. The residual interaction causes configuration 
mixing among the basis states {|/i(J))} and generates energy eigenstates \a(I)) of the 
many-body Hamiltonian, 

K')> = E*;(W)> (io) 

whose amplitudes are determined by numerical diagonalization. The diagonalization 
is done separately for each spin and parity with the lowest 700 basis states. 
The E2 operator associated with the collective rotation is given by 



(M/)|M(£2,A = 2,^ = 2)| / /(/-2)> =^^-Q 5^, (11) 



assuming that the collective E2 transition takes place between states having the same 
microscopic configuration. Static quadrupole moment Qq is assumed to be constant 
for all configurations. The E2 strength for the transition from a level a at / to a level 
(3 at I— 2 is given by 



= ^Q>^ (12) 



with normalized strength 



(13) 
(14) 




3 Level density 



In this section, we discuss characteristics of level densities of superdeformed nu- 
clei. Since our purpose is to understand mechanism of the configuration mixing and 
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Figure 2: (a)Level densities p(£ l ),(b)two-body 
level densities p2body(-E) and (c) average se- 
niorities (u) (E) of superdeformed states as 
a=A SD /20 functions of the excitation energy measured 
from the yrast line at / = 50?i in 146 Gd, and 
those for the normally deformed nucleus 168 Yb 
at I = 40h. The level density is averaged 



a=A SD /60 over both parities. The Fermi gas expressions 
are also plotted with dotted curves for var- 
ious values of level density parameter a, in 
which Asd — 150 for superdeformed states 
and A N d = 168 for 168 Yb. 



the rotational damping, we deal with the level densities before the shell model di- 
agonalization, i.e., all the calculations in this section refer to the unperturbed basis 
states and their energies {E^I)}. Fig. 0(a) shows calculated level density 

of superdeformed states and compares it with the level density in a typical normally 
deformed nucleus 168 Yb. The calculated level density of super-deformed states is sys- 
tematically much lower by an order of magnitude than that of normally deformed 
states. A similar results is obtained in the previous analysis 0, [|. This is attributed 
to the low level density of the single-particle levels at the Fermi surface (See Fig.|I|), 
which is an general feature intrinsic to the superdeformed single-particle potential. 
In order to give quantitative discussion, the calculated level density is compared in 
Fig. |2|(a) with the Fermi-gas expression [00, Q 



p(E) = Y-aT^E--* exp 2v / aE (15) 
48 

where the level density parameter a is supposed to be reduced and is treated as a 
variable parameter. The level density in the normal deformed nucleus 168 Yb is well 
fitted by the Fermi-gas expression with the level density parameter a = ^go ~ A/10, 
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Table 2: The level density parameters a to t , O2body, tt sen iority, which are obtained by a fit to the calcu- 
lated level density p(E), the two-body level density p2body(E), and the average seniority (v) (E) at 
E = 2.5 MeV and I = 50h. Since the level density parameter is often expressed as a = A/ao, the 
table lists the value of clq in unit of MeV instead of a. The fifth column lists the shell correction 

rp / 77>micro 77ismooth\ 

energy Ashcii — - — )■ 



which agrees with the standard estimate of the single-particle level density g [g(J . 
In contrast, the calculated level density of superdeformed states corresponds to much 
smaller value around a = A/15 ~ A/25. It is known from fitting experimetal data to 
the Fermi gas formula that the level density parameter in spherical closed shell nuclei 
is sizably smaller than the over- all systematicsf20| . 



Another level density function which counts the states connected by the two-body 
residual interaction plays an important role in the configuration mixing of unperturbed 
states^. This quantity, which we call the two-body level density p2body(E) in the 
following, can be calculated explicitly as described in Appendix. The result is depicted 
in Fig. ^|(b) and is also compared with the p2bod y in normally deformed 168 Yb and the 
Fermi-gas estimate for a fixed spin and parity || 

p 2b ody = |^ 5/2 £ 3/2 - (16) 

The two-body level density for superdeformed states is again much lower than that 
in the normally deformed nucleus. The corresponding value of the level density pa- 
rameter is around A/12 to A/20. 

Fig. H(c) shows that average seniority(average number of particles and holes excited 
from the reference configuration). This quantity for superdeformed states is also signif- 
icantly low in comparison with the normally deformed. The level density parameter 
^seniority which is obtained by a fit to the Fermi gas expression |21], U 

(u) = ]D S ^aE) (17) 



is very low, ranging around A/15 to A/60. The level density parameters a tot , a 2 body 
and a sen i or ity obtained from the calculated p(E), P2bod y (E), and (v) (E) are listed in 
Table. |. 

It should be emphasized that the level densities p(E) and p2body(E), and the av- 
erage seniority (u) (E) for superdeformed states are not well fitted by the Fermi-gas 




excitation energy E ex (MeV) excitation energy E ex (MeV) 



Figure 3: The level density of unperturbed np-nh configurations for (a)normally deformed nucleus 
168 Yb at / = 40h and (b)superdeformed 150 Gd at / = 50h as a function of the excitation energy 
measured from the yrast line. In 150 Gd only lp-lh and 2p-2h components are dominant for E cx < 3 
MeV, while in 168 Yb lp-lh to 4p-4h contributes to the level density at the same energy region. 

formulas with a single value of a for the super deformed states. Apparently the ex- 
citation energy dependence of the calculated level densities deviates from the Fermi 
gas formula in the excitation energy range displayed in Fig. |2|. For example, the total 
density p(E) increases very slowly for E < 1 MeV while it shows steeper increase 
at E > 1 MeV than the fitted Fermi gas formula. It is also noticed that the level 
density parameter, say atot, which fits the calculated level density p at a given exci- 
tation energy differs from the parameter embody fitting p2body and from a sen iority; «tot is 
systematically smaller than embody, and a S eniorit y being further smaller than a tot . These 
features clearly indicate that the level densities of the superdeformed states cannot 
be represented by the Fermi gas model. The deviation between a tot and a S eniorit y is 
previously pointed out[|]. 

In order to understand the deviation from the Fermi gas model, it is useful to 
classify the excited states in terms of the number of excited particles and holes. The 
decomposed level densities are shown in Fig. ^ for superdeformed states 150 Gd as well 
as for normally deformed 168 Yb. In 168 Yb, there are contributions from lp-lh to 4p-4h 
states for excitation below 3 MeV, while only lp-lh and 2p-2h states contribute in the 
same energy region in superdeformed 150 Gd. Superdeformed states with 3p-3h and 
higher np-nh configurations lie at higher energy region. Note that the shell structure 
in superdeformed nuclei produces the gap of about 1 MeV at the Fermi surface in the 
single-particle spectrum , because of which each particle-hole excitation costs at least 
about 1 MeV of excitation energy. Such a feature cannot be represented by the Fermi 
gas model which assumes uniform single-particle spectrum. 

The difference between a to t and embody is also explained by the presence of shell gap. 
The level density of np-nh states depends on excitation energy as p npn h oc E 2n ~ x for 
an equidistant single particle spectrum pl| , which is modified as p npn h oc (E — nA) 2n ~ 1 
by the presence of the gap A at the Fermi surface. The decrease of the density due 
to the gap becomes very significant for np-nh states with large n. On the other 
hand, we find from similar decomposition of the two-body level density that the major 
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contribution to p2bod y comes from the component representing the interaction of lp-lh 
to lp-lh and 2p-2h to 2p-2h states. The subcomponent associated with the interaction 
among npnh states depends only linearly on energy as p2bod y ,n P rth oc (E — nA) 1 . Thus 
the gap is less effective for p2bod y than for p tot , explaining the result shown in Fig. |2] 
and Table. that the level density parameter a tot fitting the total level density p is 
systematically smaller than a 2 bod y of the two-body density p2bod y - 

The calculated level densities of superdeformed states in A ~ 150 region are plotted 
in Fig. |3] and |^. Sizable change in the level densities depending on neutron and proton 
number N and Z is remarkable. For example, there is variation of more than factor 
2 in the total density p(E). The fitted level density parameters a to t and Gt2bod y vary 
50% as shown in Table. The variation of the level densities should be originated 
from the changes in the single-particle structure in different superdeformed nuclei, 
which are associated with the changes in the Fermi surface and the deformation. In 
order to demonstrate the correlation between the level densities and the single-particle 
structure,we compare in Table. |2|, systematics of the fitted level density parameters 
with that of the shell correction energy -E s heii — — (_E micro — ^; sm °° th ) Q f the reference 
configuration, which is a measure of the shell effects in the single-particle spectrum. 
The correlation is clear in Table. [| the larger is the shell correction energy, the smaller 
become the fitted level density parameters a to t and a 2 bod y - Taking Gd isotopes as an 
example, the level densities increase as the neutron number decreases from 86 to 82. 
Correspondingly the shell correction energy decreases for lighter isotopes. This can 
be related to both the neutron and the proton single-particle routhian structure. In 
150 Gd ,the routhian spectra are not very different from Fig. [I]. The neutron and the 
proton spectra has a large gap at N = 86 and Z = 64 Fermi surface Respectively. 
For lighter Gd isotopes ,the Fermi surface is situated at one to several orbits below 
the = 86 gap,where the single-neutron level density is rather high. Furthermore 
the proton gap at Z = 64 Fermi surface decreases as the quadrupole deformation 
decrease(from e<i = 0.556 at 150 Gd to t<i = 0.510 at 146 Gd). Both effects cause the 
increase of level density and the decrease of the shell correction energy for lighter Gd 
isotopes . Another example is 151 Dy ,in which the level density is the largest among 
the calculated nuclei. This is because there is no big gap in the neutron spectrum at 
N = 85 Fermi surface, and the proton gap at Z=66 is not very large (The proton gap 
for 151 Dy is smaller than that in Fig. [I] because of smaller quadrupole deformation). 
In a few nuclei, the calculated level densities show dependence on spin of about factor 
2. The significant spin dependence is also attributed to the changes in the shell gap. 
For example, in 151 Dy, the Z=66 proton gap just above the 64 orbit decreases with 
increasing rotational frequency up to uj ~ 0.5MeV/ft because of the signature splitting 
of 63, 64 orbits (This behavior is seen in Fig. [[]). 

4 Onset of rotational damping 

4.1 Results 

Every unperturbed many-particle and many-hole configuration created upon the 
cranked single-particle orbits would form a superdeformed rotational band if there 
were no configuration mixing. The unperturbed rotational bands are displayed in 
the upper panel of Fig. |] for the signature and parity quantum number (0, +) in 
152 Dy. The lowest band has the configuration filling the single-particle orbits up to 
the Fermi level at N=86 and Z=66 (See Fig. [I]), which corresponds to the observed 



yrast super deformed band. The second lowest band is a proton lp-lh configuration 
with the hole at the 64 orbit and the particle at [413] |. Other unperturbed bands 
displayed in the figure have mostly lp-lh and 2p-2h configurations. 

With the residual two-body interaction included, the unperturbed configurations 
are mixed to form energy eigenstates. The resultant energy levels are plotted with 
little horizontal bars in the lower panel of Fig. ^| for (0,+) states in 152 Dy. The 
configuration mixing due to the residual interaction is strongly dependent on the 
excitation energy. An eigenstate near the yrast consists mostly of a single dominant 
configuration, while a state at high excitation energy consists of many unperturbed 
configurations. Correspondingly, the E2 decay from a near-yrast level at I feeds 
to a level at I — 2 containing the same dominant component. The associated E2 
strength exhausts most of the total strength for the collective rotation. The solid 
lines connecting the levels in Fig. |] correspond to such strong E2 transitions. On 
the other hand, the E2 strength for higher excited levels is fragmented over small 
components decaying to many final states. A sequence of the strong collective E2 
transitions (solid lines) indicates a rotational band structure. The band structure 
is dominant near the yrast line while it tends to disappear as the excitation energy 
increases, indicating the onset of the rotational damping. 

In order to make a systematic analysis of the onset of rotational damping, let us 
discuss the excitation energy where the damping sets in. It is useful to define the 
branching number |6|] 

^branch («) = ^£ ( 18 ) 

associated with the E2 strength for decays from a level a. Here w a ^/3 is the normalized 
E2 strength, Eq. (|T3"D, from a level a at 7 to a level (3 at 1—2. This quantity measures 
effectively the number of E2 decay branches. The onset of the damping can be defined 
by the condition that the E2 decays have more than two branches or branch > 2 || 0]. 
An example of the calculated n branch is plotted in Fig. |7| as a function of the excitation 
energy measured from the yrast line. The branching number increases with excitation 
energy, and the onset energy of the damping E onset where n branch exceeds 2 is given as 
1.8 MeV in this example. The calculated onset energy is shown in Fig. [8] as a function 
of spin for super deformed nuclei in A ~ 150 region. 

Because of the onset of the rotational damping, there exist only limited number 
of rotational bands which show up as the regular sequences of levels connected with 
strong collective E2 transitions. The number of rotational bands thus gives a quanti- 
tative measure of the onset of the damping. Noticeably, there exist many "rotational 
bands with short length" in which strong E2 transitions continue only for two or 
three steps. In defining the number of rotational bands, we count the bands with 
length more than two. This definition corresponds to the experimental definition of 
the effective number of paths [|J which is extracted from the spectral fluctuation at 
the first ridge of the E 7l x E l2 spectra, which are formed by two consecutive E2 
transitions along rotational bands. In practice, we require that along a rotational 
band the strong collective E2 transitions continue at least two steps with the crite- 
rion rebranch < 2. Counting those levels satisfying this condition for all parities and 
signatures, say I n = 40 + , 40", 41 + , 41", we define the number of bands iV band at a 
representative spin, 40h in this example. 

The result is plotted in Fig. |[ The fluctuations in the calculated iV band depending 
on spin arise from the sharp boundary branch < 2 defining rotational bands, and 
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Figure 6: The calculated energy levels plotted with little horizontal bars for (0,+)states in 152 Dy. An 
rigid-rotor rotational energy _E r igid = 1(1 + 1)/190 (MeV) is subtracted. Solid lines connecting the 
energy levels represent strong E2 transitions with normalized strength larger than ^/l/2 = 0.707. 
The upper panel shows unperturbed levels with no residual interaction. The lower panel shows the 
results with the residual interaction. 
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Figure 7: The branching number rtbranch as a function of excitation energy above yrast line for 152 Dy, 
averaged over / = 52 ± — * 50 ± , I = 53 ± — * 51 ± . The horizontal lines shows 7ibranch = 2 used to define 
the onset of rotational damping. 

it may rather represent ambiguity of the definition of iVb an d- In order to see the 
dependence of the species,it is preferable to take averages of the onset energy E onset 
and the number of bands iVb an d over different spins. Table. |3| summarizes the average 
-^onset and iVband where the averages are taken over spin range I = 40 ~ 60h, which 
approximately correspond to spins where the superdeformed states are expected to 
dominate the yrast region. 

It should be emphasized that, although the onset energy thus defined tells approx- 
imately where the rotational damping sets in, the transition from the rotational band 
structure to the rotational damping is not very sharp as function of excitation en- 
ergy, but it develops rather gradually. As shown in Fig. there are numerous "short 
bands" which are surrounded by the excited levels with damped rotational transitions. 
Because of the presence of such scars of rotational bands the onset energy E onset 
should be regarded as one of the measures defining the transition. In fact, the scars 
still exist at 1 MeV above the onset energy. 

4.2 Discussion 

The calculated onset energy -E on set for the superdeformed nuclei is about 1.5 ~ 3.0 
MeV above yrast line, depending on species (Fig. |8] and Table. [3]). It is signifi- 
cantly higher than the onset energy -Eonset ~ 0.8 MeV in normal deformed rare-earth 
nuclei H |7J. The calculated number of bands iVband ~ 40 —70, is also larger than the 
corresponding number iVband ~ 30 in normally deformed nucleus 168 Yb, while iVband 
in a few superdeformed nuclei such as 143 Eu and 151 Dy is comparable with that in 
normal deformed nuclei. 
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Table 3: The calculated onset excitation energy -Eonsct of rotational damping measured from the 
yrast line, and the calculated number of superdeformed bands iVband, averaged over / = 40fi. ~ 60?i. 



Another characteristic feature seen from Fig. |||9] and Table. || is variation of the 
onset energy and the number of bands in different superdeformed nuclei. For example, 
the onset energy becomes the highest(about 3.0 MeV) and the number of bands is 
also large(about 70) in 150 Gd, while the same quantities become a factor 1/2 small 
in 143 Eu (about 2.0 MeV and 40) and 151 Dy (about 1.4 MeV and 40). One can also 
trace that i?onset and iVband decrease as the neutron number decreases in Gd isotopes. 

Both the difference from the normal deformed nuclei and the variation among the 
superdeformed nuclei are explained in terms of the characteristics of the level densities 
discussed in Sec. ^. In this connection, it is useful to recall the argument in Ref.0 
that the damping of rotational motion takes place if the spreading of unperturbed np- 
nh configurations is caused by the two-body residual interaction, which is represented 
by the condition T M > embody Here = 2nv 2 p2body is the spreading width of np- 
nh configurations with v being the average size of matrix elements of the residual 
interaction, and p2bod y being the two-body level density. The spreading is meaningful 
only when it is larger than the spacing between interacting states embody = PiTbody 
The above condition reads 

P2body > 1/V27TV (19) 

implying that the onset of the rotational damping is governed by the ratio between 
the two-body level spacing embody = P^body anc ^ ^he two-body interaction strength. 
As P2bod y (-E0 increases with excitation energy, the rotational damping sets in at the 
excitation energy where the condition Eq. (19) is satisfied. It is now to be noticed that 
the two-body level density in superdeformed nuclei is much lower than in normally 
deformed nuclei as shown in Fig. 0. Assuming that the matrix elements of residual two- 
body interaction does not depend on deformation, the onset energy E onset satisfying 
the condition Eq. fll9] ) must be much larger in superdeformed nuclei than in normally 
nuclei. The variation of .Sunset among the superdeformed nuclei can be explained 
along the same argument. It is seen from combination of Table. ^| and Table. |^ that 
the variation of -E ori set is correlated with that of the fitted level density parameters 
and the shell correction energy. For example, the onset energy E onset , the fitted level 
density parameters and the shell correction energy become the largest in 150 Gd. In 



this nucleus, the two-body level density becomes very small due to the shell gap at 
N=86 and Z=64 Fermi surfaces. 

To be more quantitative, let us adopt the Fermi gas expression Eq. ( |i~6|) for the 
two-body level density. Then the onset energy is estimated from Eq. ([19]) as 



£ onset = 7.1KT 2 / 3 a- 5 / 3 . (20) 

Choosing v = 25 keV and a = A/10(A = 170), the onset energy in this estimate 
E nset — 0.74 MeV agrees with the microscopic calculation for normally deformed 
nuclei 0. With the same interaction strength, but with smaller level density parameter 
a = A/ 15 ~ A/20(A = 150) corresponding to the superdeformed states, the estimate 
reads -E nset — 1-8 ~ 2.9 MeV, which roughly agrees with the microscopic results in 
Table. 

Let us now discuss the number of bands. The microscopic results are about 40 ~ 70, 
which are about factor 1-2 larger than in normally deformed nuclei. The variation of 
iVband among the superdeformed nuclei is correlated to that of the level densities as 
seen from Table. [| and Table. § the larger the number of bands is, the smaller are the 
level density parameters. For qualitative understanding of these features, let us utilize 
two schematic models. The first one is along the same line with the above discussion 
for i? onS et- The number of bands in this model may be estimated by counting the levels 
up to the onset energy given by Eq. (p0|) provided that the boundary is sharp(while,in 
fact, this assumption is not necessarily correct as Fig. || indicates). Then the number 
of bands is given by 



iV band = 4 / p{E)dE 
Jo 



! ^ J X ° aT 5/4 exp 2^/x~dx, (21) 
x = 7.11(cw;)- 2/3 . 



Here the Fermi gas expression Eq. (|I5D for the total level density p(E) is used. Multi- 
plier 4 is needed to count different parity and signature states. The number of bands 
is an steeply increasing function of (av)^ 1 ; 

AWd = 26, 55, 100, 170, 270, 400 for (av)~ l = 2, 3, 4, 5, 6, 7. The choice of the interac- 
tion strength v = 25 keV gives iVband — 34 for a normally deformed nucleus with A = 170 
in agreement with the microscopic calculation. With small level density parameters 
corresponding to the superdeformed nuclei, the estimate produces iV band =101 ~201 
for a — A/15 ~ A/20. This indicates the significant increase of the number of su- 
perdeformed bands and its sensitivity to the variation of level density. However this 
model apparently overestimates the microscopic results in Table. || The overestimate 
in iV band can be traced back to the fact already revealed that the density parameter a 
fitting the total level density p(E) is significantly smaller than that for p2body{E). As 
is discussed in Sec. |3|, the level densities affected by the presence of the shell gap in 
the single particle spectrum cannot be simulated simply by a single value of the the 
level density parameter. In order to treat the shell gap explicitly, let us next consider 
the other model, which simplifies the situation by assuming that the single-particle 
spectrum is uniform with density \g Q (related to the level density parameter a— \go) 
except the presence of the gaps with interval A at the Fermi surface. Considering the 
limit of large A, the lowest excited levels are dominated only by lp-lh excitations. In 



this limit, the total level density p(E) is significantly reduced, and is expressed as 



p(E) = P 2 body (E) = l -gl{E -A) (E > A) (22) 

The two-body level density is the same as p(E) since all the lp-lh states can interact 
through the two-body force. With v = 25 keV and a = A/10 corresponding to the 
average single-particle density pD|, the onset energy and the number of bands are 
estimated as 

4\/2 

£ onset = —= — - + A = 1.53MeV + A, (23) 

iVband = ^-^= 49. (24) 

The onset energy E onset increases with increase of the single-particle gap A while the 
estimated number of bands reaches to a constant independent of A, which is not very 
large comparing with an estimate iVband = 44 corresponding to A = 0(a = A/10, A = 
150 in Eq. (0))- The two schematic models discussed above simplify the realistic 
features of the single-particle spectra for superdeformed potential. The calculated 
level densities presented in Sec. || seem to lie in-between the two models. 

It should be reminded that the onset of rotational damping does not take place 
very sharply at the onset energy, but proceeds quite irregularly as function of the 
excitation energy, as is discussed in the preceding subsection and demonstrated in 
Fig. H This is because there are local fluctuation of the unperturbed levels and 
the matrix elements of the residual interactions depending on the configurations of 
the many-particle many-hole states. A little deviation of -E on set and iV band from the 
systematics expected from the N, Z-dependence of the shell structure effect in the 
single-particle levels may be due to such fluctuations. 

We noticed, in Fig. |9], systematic spin dependence that the number of bands in- 
creases with decreasing spin below 40h. This spin dependence is not correlated with 
the level density since the calculated level densities in Fig. |4] and [5| have no systematic 
spin dependence which can explain the increase of iVband- This may indicate that the 
onset of rotational damping is not controlled by the configuration mixing alone. It 
is should be noticed in this connection that the alignment effect is also important 
to cause the rotational damping, while it is not very large for the low spin region in 
superdeformed nuclei. The dispersion of rotational frequency 2Au;,which measures 
the alignment effect and represents how different rotational bands respond to the spin 
change of A J = 2, is about 40 keV at / = 30 on the basis of the estimate of Ref. 0. 
This is not sufficiently large than the scale of the residual interaction effects ~ 50 
keV (embody ~ V2nv at the onset energy), the rotational correlation persist more than 
expected from the argument based on the level densities. It is to be reminded that pos- 
sible persistence of the rotational band structure in the region of strong configuration 



mixing is pointed out in Ref. [23 



4.3 Comparison with experiments 

Recent experiments identify several superdeformed bands in a nucleus in A ~ 150 
region [p]]. Our microscopic calculation predicts at least more than ten superdeformed 
rotational bands even if "the short rotational bands" are excluded (there are three 
other figures in addition to the lower panel of Fig. ^ for 152 Dy). It implies that 




Figure 10: The experimental effective number of decay paths [|l(J and the calculated number of 
superdefomed bands N^and as a function of the transition gamma-ray energy. 



the current spectroscopic studies of discrete gamma-ray peaks reaches far below the 
onset of rotational damping in superdeformed nuclei. In other words, the present 
calculation suggests possibility to find further number of superdeformed bands in a 
nucleus in future experiments. 

Available experimental data concerning the the onset of rotational damping in su- 
perdeformed nuclei come from analysis of quasi-continuum part of double-coincident 
E lx x E 12 spectra although the number of data is very limited so far. An indirect 
information comes from an simulation of the quasi-continuum gamma-rays in superde- 
formed 152 Dy pi| , which points to an estimated value -E on set around 2.5 ~ 3.0 MeV 
from the observed intensity of the ridge part of the i? 71 x E 12 spectra [pSfl . This anal- 
ysis agrees with the present microscopic calculation in the qualitative feature that 
the onset of damping takes place at much higher excitation energy than in normal 
deformed nuclei. 

More direct information comes from the fluctuation analysis method which extracts 
the effective number of paths of decays from the quasi-continuum part of E lx x E l2 
spectra [§, £|. When this method is applied to the first ridge of the spectra located 
at £7 7l — E l2 = ±4/X lines with X being the moment of inertia of superdeformed 
bands, the effective number of paths counts the number of two consecutive collective 
E2 transitions, 1+2 — >/— > 1—2. Thus the effective number of paths at the first ridge 
approximately corresponds to the number of bands which is defined in Sec. |2|. The 
data for superdeformed nuclei is available only in 143 Eu Jl(J. Comparison between the 



theory and the experimental data is done in Fig. [10], in which our results is plotted as 
a function of the average gamma-ray energy. The data points are limited in the region 
E 1 ~ 1250 — 1550 keV where the quasi-continuum ridge is observed. The calculated 



number of bands agrees with the experimental effective number of paths at around 
£ 7 ~ 1400keV. 

The effective number of paths becomes smaller around high and low ends of the 
measured E 1 range than the theoretical number of bands. This may suggest some 
other mechanisms that destroy rotational band structures in addition to the rota- 
tional damping. For example, for the region with <; 1500 keV corresponding to the 
very high spins / <; 55/i, the fission decay may compete with the collective E2 gamma 
decay and hinder the E2 decay branches^]. Furthermore , feeding mechanism of the 
superdeformed states may affect the effective number of paths. The low gamma-ray 
energy E 1 ^ 1300 keV, on the other hand, corresponds to spins I ^40h, at which the 
lowest superdeformed states lie higher in energy than the yrast line consisting of nor- 
mally deformed states. Since the superdeformed states are surrounded by numerous 
normal deformed states in this situation, mixing between normal and superdeformed 
states may take place and cause breaking of rotational band structures of the superde- 
formed states. Inclusion of these effects has possibility to explain the data points, but 
it is beyond the scope of the present paper. 

5 Conclusion 

We analyzed the onset of the rotational damping in A ~ 150 superdeformed nuclei 
by means of the shell model diagonalization method comprised of the cranked Nilsson 
mean-field and the residual two-body interaction. 

The rotational damping in superdeformed nuclei sets in at the excitation energy of 
about 2 ~ 3 MeV above superdeformed yrast line. The onset energy is much higher 
than the theoretical prediction in normally deformed rare-earth nuclei. The effective 
number of rotational bands is calculated to be 40 ~ 70, which is slightly larger than 
the corresponding number (about 30) in normally deformed nuclei. Both the onset 
energy and the number of bands show the variation depending on the superdeformed 
species even among A ~ 150 region. 

All these features characteristic for the superdefomed nuclei are originated from 
the shell structure in the single-particle spectrum associated with the superdeformed 
shape of the mean-field potential. This was demonstrated explicitly by looking into 
the total level density as well as the two-body level density which is relevant to 
the configuration mixing via two-body forces. The shell gap in the single-particle 
spectrum causes significant decrease in the level densities in comparison with those 
in the normally deformed nuclei, and causes increase in the onset energy and the 
number of bands. Since the shell structure at the Fermi surface varies with neutron 
and proton numbers and associated changes in deformation, the level densities show 
significant dependence on superdeformed species even among A~ 150 region. It causes 
the variation of the onset of rotational damping. 

The calculated results are consistent with existing data. More systematic com- 
parison with data is required, however, in order to check the predictions such as the 
mass-number dependence, and also to clarify other mechanisms which are not included 
in the present theory. 
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Appendix 

In order to calculate the two-body level density p 2 bod y (E), it is useful to define a 
related quantity 

( El , E 2 ) = £ '5{E 1 - E^)5(E 2 - Eft (25) 
mm' 

where Yluu' ' denotes the summation over pairs of np-nh states which interact via 
two-body residual interaction. Eq. (p5|) can be rewritten as 

pQ\E u Ek) = £<*(£!- £ n ) E 'S(E 2 -E I1/ ) (26) 

M M '0) 

= E 5 (^l-^)P2body(/i,^ 2 ) (27) 
M 

= p(E 1 )p 2hody (E 1 ,E 2 ) (28) 

Here the quantity p 2 body(p, E 2 ) representing density of states coupled with the state 
p via the two-body interaction is assumed to depend only on the energy of the 
state. Then the two-body level density is given by 

P2boAy{E) = p 2 body(E, E) 

= pW(E,E)/p(E) . (29) 

In numerical calculation, the delta functions in Eq. (]25f) is replaced by an Gaussian 
function used in the Strutinsky averaging. The level density p(E) = 8(E — E^) is 
also calculated with the same Gaussian function. The width of the Gaussian function 
used in the numerical calculation is 200 keV. 
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